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ABSTRACT 

We report the discovery of 6.7 GHz methanol maser periodic flares in four massive star form¬ 
ing regions and the updated light curve for the known periodic source G22.357-1-0.066. The 
observations were carried out with the Torun 32 m radio telescope between June 2009 and 
April 2014. Flux density variations with period of 120 to 245 d were detected for some or all 
spectral features. A variability pattern with a fast rise and relatively slow fall on time-scale 
of 30-60 d dominated. A reverse pattern was observed for some features of G22.357-1-0.066, 
while sinusoidal-like variations were detected in G25.411 -hO. 105. A weak burst lasting ~520 d 
with the velocity drift of 0.24 km s ' yr * occurred in G22.357-1-0.066. For three sources for 
which high resolution maps are available, we found that the features with periodic behaviour 
are separated by more than 500 au from those without any periodicity. This suggests that the 
maser flares are not triggered by large-scale homogeneous variations in either the background 
seed photon flux or the luminosity of the exciting source and a mechanism which is able to 
produce local changes in the pumping conditions is required. 
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1 INTRODUCTION 


The 6.7 GHz methanol maser transition has been detected from 
very many massive star-forming regions JPandmnetaljj2007|; 


very many massive star-tormmg regions JPanomnetalJJZtHM 

Caswell et al.l201(l[20ni:lGreen et al.l2O10Ll2O12l:ISzvmczak et "sTl 


2012[) . Since the early report of icaswell et alJ h995ll . the tempo- 

ral properties of methanol masers have been intensively studied. 
In particular, the phenomenon of periodic or regular variability at 
6.7 GHz in high-mass young stellar objects (HMYSOs) is a matter 
of active debate. 

To date, periodic variations in the maser emis- 
sion have been detected f r om a do z en HMYSOs 
([Goedhart. Gaylard & van d er W ald | 2001^ |200 4 iGoe dhart et alJ 
2d09[ l26T3r~TArava et alflT^ ISzvmczak et al.l |201ll: 


Fujisawa et alJ 12014 1. They are characterized either by flaring 


or sinusoidal patterns with periods ranging from 29 .5 to 668 d and 
the re lative amplitude from 0.2 to more than 30 jFuiisawa et al.l 
|2014) . Some of the sources show periodic behaviour in part of 
the spectrum, whereas others display synchronised variations 
of all maser feature s, sometimes with de l ays in flaring betwee n 
individual features iSzvmczak et al.l[20lll : IGoedhart et al]|2014h . 
The 12.2 GHz maser emission, when observed, commonly shows 
periodic amplitude variations higher than those at 6.7 GHz 
(IGoedhart etS]|20l4 . G9.62+0.20E, an archetype of the class, 
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exhibits a re gular flaring behav i our al so in the 107 GHz methanol 
maser line ( Ivan der Walt et al] 1200911 . Correlated periodic vari¬ 
ability between the 6.7 GHz methanol and 4.8 GHz formaldehyd e 
maser lines has been reported for G37.55-0.20 l lArava et al.l201(]h . 

Different models have been proposed to explain the regular¬ 
ity of the methanol maser flares. For sources with periods of 240- 
400 d it has been postulated that a binary interaction between a 
HMYSO and its companion could periodically alter either the seed 
photon flux or pump rate in the maser region ivan der Walt et all 
l2009l : Ivan der Waltiboilh . Specifically, a colliding-wind binary 
(CWB) could produce a periodic pulse in the ionising radiation 
that changes the electron density and periodically modulates the 
free-free emission from the Hii region against which the maser 
clouds are projected. In another model, the periodic maser varia¬ 
tion is related to the flux of infrared emission that could be reg¬ 
ulated by periodic accretion of material from the c i rcumb inary 
disc o nto a protostar or accretio n disc Arava et alj l201(]|l . Re¬ 
cently, IParfenov & Sobole^ l l2014ll proposed that the dust temper¬ 
ature variations in the circumbinary accretion disc are caused by 
rotation of hot and dense material of the spiral shock wave in the 
disc central gap. The maser brightness depends on the density of 
the shocked material along the line between the central star and 
the maser region. Simulations of stellar evolution indicate that a 
HMYSO still accreting at a rate larger than lO^^MQyr^’ and with 
a radius exceeding l OORfri at maximum is pulsationally unstable 
dlnavoshi et alJl2Rl4 . The expected variations of luminosity with 
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typical periods of 10-100 d could drive the periodic variability of 
met hanol masers wi t h a sin usoidal light curve. 

lOoedhart et alJ j2005h mapped G9.62+0.20E at 12.2 GHz us¬ 
ing Very Long Baseline Interferometry (VLBI) and found that the 
brightness of the maser components increased during the course 
of a flare, whereas the overall structure was unchanged. The mor¬ 
phology of the 6.7 GHz maser emission of the periodic source 
G22.357-1-0.066 was also unchanged in both observations near the 
burst maximum dSzvmczak et al.ll201 ih . These findings imply that 
the causes of the flares originate outside the maser region. 

In this paper we report the discovery of periodic variations in 
four 6.7 GHz methanol sources and produce an up dated light curve 
of th e known periodic source G22.357-1-0.066 l lSzvmczak et al.l 
l201lh obtained as part of a long-term monitoring programme of 
a large sample of obje cts detected with the Torun 32 m telescope 
dSzvmczak et alJboi^ . 


2 OBSERVATIONS AND DATA ANALYSIS 

134 target sources were selected from the Torun methanol maser 
catalogue dSzvmczak et al.l2012h using coarse criteria intended pri¬ 
marily to monitor the objects with peak flux density greater than 
5 Jy and declination north of -11°. These objects stay above the 
elevation limit for relatively long durations and are suitable for a 
survey with the Torun 32 m radio telescope. A detailed description 
of the sample will be presented in a subsequent paper. 

Each object was observed at least once a month during the 
period from 2009 lune to 2010 October; then the objects showing 
a strong variability were observed at irregular intervals of 5-lOd, 
until 2014 April. Several gaps of 3-4 weeks in the monitoring cam¬ 
paign occurred due to scheduling constraints. 

The 6668.519 MHz methanol maser transition was observed 
in a frequency-switching mode using a dual-channel HEMT re¬ 
ceiver mounted on the 32 m radio teles cope. The detail s of th e 
observing setup have been described by ISzvmczak et al.l ( l2014h . 
In summary, the system equivalent flux density was about 300 Jy, 
the half-power beam width 5.5 arcmin and the pointing accuracy 
22 arcsec. The autocorrelator was configured to record 4096 chan¬ 
nels over 4 MHz for each polarisation, yielding a spectral reso¬ 
lution of 0.09 km s“* after Hanning smoothing. The typical rms 
noise level in the spectra was 0.20-0.25 Jy after averaging both 
polarisations. The system was calibrated b y measuring the con¬ 
tinuum point source 3C123 dOtt et al.lll994h . Further checks were 
undertaken using the methanol maser source G32.745-0.076, in 
which some spectral features have not shown variability higher than 
~8 per cent dSzvmczak et alj2014h . The errors on the absolute flux 
density measurements, estimated from the standard deviation of the 
flux density values of ’non-variable’ features of this source, were 
usually less than 10 per cent. 

The edited and flagged light curves were typically de-trended 
using a third-order polynomial or lower prior to period search. Fig¬ 
ure [T] shows the effect of these procedures for the light curve of 
80.09 km s“* feature of source G22.357 -f 0.066. Data points were 
discarded if their flux densities differed by ^ 5<t from the values 
observed just on the previous and the next day. This criterion was 
applied only for the intervals of quiescent emission. In the flare 
intervals sampled daily the individual observations were discarded 
only if the flux density obviously dropped below ^ 5cr as compared 
to the adjacent points. In effect about 2 per cent of observations of 
G22.357-)-0.066 were discarded while for the rest sources the per¬ 
centage of flagged observations is smaller. Slow and low amplitude 
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Figure 1. Effect of flagging and de-trending of the data on an example light 
curve. The light curve of 80.09 km feature of source G22.357 -f 0.066 is 
shown before flagging (a), after flagging but before de-trending (b) and after 
flagging and de-trending (c). The dotted horizontal line in panels a) and c) 
is put to mark the difference in the quiescent levels of the original and final 
curves. 

changes of the flux density in the quiescent state were fitted with a 
third-order polynomial (Fig. [Til. In several cases a fit with a linear 
function was enough. 

In order to calc ulate the per i od and assess its significance 
the LS periodogram jLom 9 ll976l;IScargi3ll982h was used. This 
well established method is effective for unevenly spaced data as 
ours. For each periodically variable feature we obtained a nor¬ 
malised phase-folded flare profile by dividing the time-series into 
cycles based on its best-fit period and the flux density values in 
each cycle were divided by the maximum observed in the cycle. 
We fitted each normalized light curve with a power function of 
the form; 5 (t) = A®*'* -F C, where A and C are constants and 
s{t) = [bcos(a)t -F — fsin(a>t -F (j)))\ -F a. Here, b is the am¬ 

plitude measured relative to the average value a, co = 2nlP, where 
P the period, (p the phase, and / the asymmetry parameter defined 
as the rise tim e from the minimum to the maximum flux, divided 
by the period dSzvmczak et al.ll201 ih . Time delays in flaring be- 
twee n features were calculat ed with the discrete correlation func¬ 
tion tedelson & Krolikl[l988h . 


3 RESULTS 

We have detected periodicity in five sources of which four are re¬ 
ported for the first time. Their average 6.7 GHz spectra are shown 
in Figure 12 and the plots of the flux density as a function of time 
and velocity are presented in Appendix A (Figs. lAT4A5t . The light 
curves of features with periodic behaviour are shown in Figures 
12 and [4] and their parameters are given in Table [T] Because of the 
non-Gaussian nature of many of the spectral features their flux den¬ 
sities were read off from the spectra. We list the velocity of the 
feature, the mean flux density, the relative amplitude, the period, 
the time-scale of variability, the ratio of the rise time to the decay 
time (Rrd), and the time delay of the flare between individual fea¬ 
tures. The relative amplitude is equal to (Smax - Sminj/Smin, where 
Smax and S,nin are maximum and minimum of the flux density, re¬ 
spectively. The time-scale of variability is the full width at half¬ 
maximum (FWHM) of the flare. 
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Figure 2. Spectra of 6.7 GHz periodic methanol masers obtained at the po¬ 
sitions given in Table[T] The solid (black) lines show the average spectra for 
the whole monitoring observations, while the dashed (blue) and dotted (red) 
lines show the average spectra for intervals of high and low emission levels, 
respectively. The flaring features are marked by vertical bars. The velocity 
intervals of the emission from nearby but significantly (<78”) offset sources 
are denoted by brackets. 


3.1 G22.357-h0.066 

In the velocity range from 78.95 to 83.5 km the emission 
showed periodic variations. Five permanent features with a mean 
flux density of 0.8-12 Jy and relative amplitude of 2-7 varied with 
a period of 177.4-179.3 d (Fig.[3] Table[T). The mean period was 
178.2±1.9d, where the error is a measure of the period dispersion 
for the features. A harmonic series is seen in the periodograms for 
all five features (Fig. lA6t . The normalised phase-folded flare pro¬ 
files are shown in Figurej^ The ranges of their FWFIM and are 
28.6-32.7 d and 0.68-1.83, respectively. The delays between the 
flare peaks of these features were 1-9 d (Table[T). 

Some features showed long-term changes of flux density in the 
quiescent state. For instance the 78.95 km s“* feature initially was 
below the level of significant emission and then increased to 0.9 Jy, 
whereas the feature 79.52 km s“' decreased from 16 to 0.8 Jy over 
the monitoring period. For these features a systematic increase or 
decrease in the peak flux density occurred in consecutive cycles. 
In other periodic features the peak flux density commonly changed 
irregularly from one cycle to another. 

The flux density of feature 83.42 km s“* in the quiescent state 
dropped below a typical 3(t sensitivity level of 0.6 Jy. We have 
made no attempt to estimate the intensity in this state with long 
integration observations, thus the relative amplitude (Table[T) gives 
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Figure 3. Light curves of selected 6.7 GHz methanol maser features in 
G22.3574-0.066. The dotted line in each panel marks the 3(T level. 


a lower limit. Faint emission with an intensity of 0.9-1.5 Jy at max¬ 
ima appeared after MJD 55885 at 78.64 km s“' and marginally 
detected emission (2.5-3.5 (t) near 82.60 and 86.50kms“' also 
showed flares synchronised with the flares of the main features 
(Fig.iMJ. No periodic variations were seen for the features near 
85.00 and 88.43 km s^*. 

After MJD 56085 new emission appeared near 85.75 km s“* 
reaching a maximum of ~1.3Jy near MJD 56248 (Fig. [3. This 
burst lasted about 520 d with no evidence of flux density modu¬ 
lation synchronised with the flaring activity of the permanent fea¬ 
tures. The line-width of the feature was 0.26-0.32 km s“* and un¬ 
correlated with the flux density. The peak velocity showed a drift 
of 0.24 km s“* yr“‘ (Fig. The event character istics are similar to 
those reported in Cep A jSzvmczak et al.ll201^ but the signal-to- 
noise ratio in the spectra is less than 6 and it is not clear whether 
the velocity drift is ca used by variability of a close pair of spec- 
trall y blended features dSzvmczak et al.ll2014l) or by motion of the 
gas iooddi et al.l201 ih . 


3.2 G25.411-h0.105 

Periodic changes in the flux density of both spectral features 
occurred until MJD 56208, then the red-shifted emission at 
97.10 km s“* slightly decreased without any periodic modulation, 
whereas periodic changes of the blue-shifted emission centred at 
94.78 km s“* were seen throughout our monitoring (Figs. |4] and 
IA2l l. In the MJD interval of 55013-56208 the period was 245 d (Ta¬ 
ble [Til. The light curve of each feature has a sinusoidal-like shape 
with Rj;] of 0.9-1.0 but the FWHM of the flare at 94.68 km s“' is 
approximately double that of the 97.10kms“* (Fig. [7] Table[T). 


3.3 G45.473+0.134 

The velocity range of the maser emission in this 
source is 11.2 km s“' (Fig. |2). The high resolution 
maps revealed that the emission in the velocity range 
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Figure 4. Light curves of flaring 6.7 GHz methanol maser features in 
sources a) G75.76+0.34, b) G73.06+1.80, c) G45.473+0.134 and d) 
G25.411+0.105. The dotted line in each panel marks the 3cr level. 



55.80-58.28 km s“’ comes from the source G45.493+0.126 
jBartkiewicz. Szymczak & van Langeveldel 12014|) that is within 
the 32m antenna beam when pointing towards G45.473+0.134. A 
clear periodicity of 195.7 d is seen for the feature with the peak 
at 59.60 km s * (Figs. |41 lA^ and lA6l Table [T). The light curve of 
this feature is characterised by irregular changes of the peak flare 
from 3.6 to 8.6 Jy, and a very stable flux density of 1.3-1.4 Jy in 
the quiescent state. The FWHM is 34 d and Rrd=0.46 (Fig.|7] Table 
[T]l. The features in the velocity range from 64.2 to 66.3 km s“* 
show weak flares synchronised with the flares of the 59.60 km s“* 
feature but only in some cycles. The emission at 62.23 km s“* is 
constant within less than 10 per cent over the whole monitoring 
period. 


3.4 G73.06+1.80 

Features at -2.59 and 5.93 km s“* showed variations with a pe¬ 
riod of 160 d (Figs. l4l |A4I and lA6l Table [T). The weak emission 
at 6.95 km s“' flared synchronously with the strongest features in 
some cycles only. The flare peak flux density varied significantly 
from cycle to cycle for all features. In the quiescent state the emis¬ 
sion of feature -2.59 km s“' increased from 0.5 to 6Jy, whereas 
that of feature 5.93 km s“' was constant within 3cr. The light curves 
of the periodic features were strongly asymmetric (Fig.|7) with R,.d 
of 0.27-0.35 and the time-scale of variability of 60-66 d (TablellJ. 
The mean delay between the flaring features was 13.8 ± 3.6 d. 


3.5 G75.76+0.34 

The maser emission at velocities lower than -9kms“' exhib¬ 
ited a flaring behaviour with a period of 119.9 d. The intensity of 
flare peaks was clearly modulated at least over the first 12 cycles 
since MJD 55158 and the periodograms for features -10.93 and 
-9.44km s“* revealed less significant (p<0.005) periods of 348.8 
and 328.6 d, respectively (Figs. I1|A5] and IA6t . The flares were 
asymmetric with Rrd of 0.63-0.72 and their FWHM were about 
33 d (Fig. a Table [T]l. The emission of feature -9.44 km s *, in 
some quiescent intervals, dropped below a typical 3 <t sensitivity 
level of 0.7 Jy, so that the relative amplitude of 7.6 (Table[T) gives 
the lower limit. The flux density of feature -10.93 km s“' always 
dropped below this threshold. No attempt was made to determine a 
low state flux density with long integration observations. The emis¬ 
sion in the velocity range from -3.2 to 1.0 km s“' did not show pe¬ 
riodic variability. The JVLA dat^i] imply that this emission comes 
from the source G75.782+0.342 which is located ~76” north-east 
from the maser source G75.761+0.340 showing periodic variations 
of both features -10.93 and -9.44kms“‘. 


4 DISCUSSION 
4.1 Period statistics 

A sample of 134 methanol masers with a peak flux density greater 
than 5Jy and declination n orth of-11° selected from the Torun 
methanol maser catalogue dSzvmczak et al.ll2012h was monitored 
with the 32 m radio telescope. We found that about 4 per cent 
(5/134) of the sources show periodic variations with periods rang¬ 
ing from 120 to 245 d. In order to examine the sensitivity of our 
observations to periodic sources we have injected periodic sig¬ 
nals of different levels and periods into real light curves and es¬ 
timated the reliability with which they are appeared in the peri¬ 
odograms. For the light curve of G22.357+0.066 consisting 279 
observations a significant frequency (a false alarm probability of 
0.01 per cent) in the periodogram occurred for signals with the 
relative amplitude >0.26 and period from 20 to 500 d. In the case 
of source G45.473+0.134 with only 86 observations a periodic in¬ 
jected signal of 30-500 d was significant for the relative amplitude 
>0.41. We conclude that our observations with average cadences of 
~6d (G22.357+0.066) and ~21 d (G45.473+0.134) are complete, 
in terms of ability of periodic maser detection, for the relative am¬ 
plitudes higher than 0.26 and 0.41, respectively. Our detection rate 
of periodic masers is a factor of three smaller than those reported by 
iGoedhart et al.l ( l2004b who found 7 periodic sources from their 54 
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Table 1. The parameters of periodic maser features. 


Source 

name 

RA(2000) 

(h m s) 

Dec(2000) 

(O , 

Feature 

velocity 

(kms“‘) 

Mean 

flux density 

(Jy) 

Relative 

amplitude 

Period 

(day) 

Time-scale 
of vaiiability 
FWHM (day) 

R.d 

Time 

delay 

(d) 

G22.357+0.066 

18 31 44.1 

-09 22 12.3 

78.95 

0.92(0.04) 

3.5(0.4) 

179.3 

32.7(1.7) 

1.83(0.39) 

0.85(0.55) 




79.52 

4.56(0.20) 

4.2(0.6) 

177.4 

28.9(1.1) 

1.50(0.28) 

1.81(1.14) 




80.09 

11.19(0.22) 

2.0(0.1) 

179.3 

29.8(1.0) 

0.68(0.08) 

0 




81.45 

2.60(0.12) 

7.0(0.7) 

177.4 

32.1(1.5) 

0.76(0.12) 

-7.00(3.38) 




83.42 

0.78(0.03) 

>4.2(0.3) 

177.4 

28.6(1.6) 

1.08(0.27) 

-4.75(2.59) 

G25.411-f 0.105^ 

18 37 16.9 

-06 38 30.5 

94.68 

7.34(0.20) 

1.2(0.1) 

244.7 

121.7(2.1) 

1.00(0.18) 

- 




97.10 

11.19(0.50) 

2.1(0.!) 

244.7 

68.9(1.8) 

0.91(0.13) 

- 

G45.473-f0.134 

19 14 07.3 

11 12 16.4 

59.60 

2.44(0.18) 

3.5(0.6) 

195.6 

34.0(1.2) 

0.46(0.10) 

- 

G73.06-fl.80 

20 08 10.1 

35 59 24.7 

-2.59 

4.67(0.18) 

1.9(0.6) 

159.0 

59.5(2.4) 

0.27(0.18) 

13.8(3.6)'’ 




5.93 

1.45(0.05) 

2.3(0.3) 

159.0 

66.1(2.0) 

0.35(0.15) 

0 

G75.76+0.34 

20 21 40.1 

37 25 37.0 

-10.93 

0.52(0.03) 

>5.4(0.8) 

119.9 

33.1(1.9) 

0.63(0.14) 

- 




-9.44 

2.00(0.14) 

>7.6(1.3) 

119.9 

32.5(1.3) 

0.72(0.11) 

- 


“ the parameters derived for MJD interval of 55013—56208, ^ the mean for cycles 1-3, 6 and 11 
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Figure 5. Phase-folded normalized light curves of selected features of 
G22.356-1-0.066. Each panel is labelled with the peak velocity. 


targets. This significant difference is likely due to a selection effect; 
Goedhart et al.’s sample basically relied on a list of targets previ¬ 
ously known as variables, while we studied a flux limited sample 
without prior information on the variability. 

Figure [8] shows the period distribution for t he 16 periodic 
metha nol sources known to date as compiled by iFuiisawa et al.l 
(I 2 OI 4 I 1 (th eir table 2) supplemen ted by our new results and refined 
data from iGoedhart et ^ ( l2014h . The period ranges from 29.5 to 
668 d with a median value of 210 d. Our detections have a rather 
narrow range of period from 120 to 245 d. The period distribution 
may suggest that the intrinsic fraction of periodic masers decreases 
with the period length. 


4.2 Maser flare profile 

We note that 9 out of 16 periodic masers are objects where 
the representative feature has an asymmetric flare profile. 


i.e. with a fast rise and relatively slow fal l of intensity. 
These are th e following sources G9.62-l-0.20E ( IGoedhart et al.l 
I200I |2014 , G22.357-f0.066, G45.473-I-0 134, G73.06-fl.80 , 
G75.76-f0.34 (Table [p, Gp.55 -f00.20 dArava et alj l2010h . 
G107 29-f5.64 dFuiisawa et al.ll20l4 . G328.24-0.55 and G331.13- 
0.24 dGoedhart et al.ll20l4 . The asymmetric flare profile can b e 
readily explained with the simple CWB model dvan der Walll201 ih . 
In this model a periodic pulse of ionizing radiation enhances the 
electron density in a volume of partially ionized gas against which 
the maser clouds are projected. The intensity of maser emission 
varies synchronously with the periodic variations of the back¬ 
ground radiation amplified by the maser process. The different 
shapes and amplitudes of the maser flares in individual sources can 
be satisfactorily re produced by modifying the orbital param eters of 


be satistactorily re produced by modity i 
the binary system dvan der Walfcoi ll ;t 


Goedhart et af 


f aram eti 

2 oI 4 . 


The mean profiles of the flares in G45.473-f0.134, 
G73.06-fl.80 and G75.76-f0.34 (Fig. |7]l are generally consis¬ 
tent with the CWB model. The changes in flare profiles between 
cycles and features of the two last objects (Fig.|4j may suggest that 
the seed photon flux is not constant and the maser optical depth 
shows long-term changes. The variability patterns of the features 
at 80.09 and 81.45kms-‘ with R,d ^0.72 of G22.357-f0.066 
are s imilar to those reported for G9.62-f0.20E dGoedhart et al.l 
I2OO9I) but they are quite different from those of the features 78.95 
and 79.52 km s“* with Rrd >1.50 (Eig. IS}. Such widely different 
patterns of the features cannot be explained by changes in the 
backgr ound free-free emi ssion postulated in the simple CWB 
model dvan der Wallll20llh . An alternative explanation proposes 
a mo del of the circumbinary accretion disc dParfenov & Sobole^ 
|2014 in which the pattern of the maser flare follows that of the 
gas column density along the line between the central star and the 
mase r region. Eor the specific structure of the disc dSvtov et al.l 
I2OO9I) . the column density changes with the inclination angle of the 
binary system, so the flare profile with R^d <1.0 may be generated 
in the molecular disc seen edge-on (i=90°). A small change of 
disc inclination on the line of sight (i=83°) would result in a flare 
profile with R^d >1.0. This issue needs detailed calculations and is 
beyond the scope of the paper. 


G25.411-f0.105 has significantly different light curves for the 
two main features (Eig.|4). The maximum of the 94.68 km s“* fea¬ 
ture leads in phase by ~ 0.15P relative that of 97.10 km s“' feature 
while both features have a minimum at the same phase. Similar 
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MJD 

Figure 6. Maser burst in G22.357+0.066. The contours are at 2, 3, 4, and 
5 times the Icr noise level of 0.26Jy. The vertical bars in the bottom axis 
correspond to the dates of the observed spectra. 

characteristics were reported for G12.89+0.49 and G331.13-0.24 
dGoedhart et al.l2014l) and could be the result of an amplification of 
the background emission from a processing outflow or thermal jet 
both of which are not directly asso ci ated with the metha nol maser 
clouds dMacLeod & Gavlamll996h . iHofner et~^ ( 1201 ih detected 
in G25.411+0.105 a weak compact continuum source at 25.5 GHz 
with a spectral index of 0.6 between 19.3 and 25.5 GHz at the west¬ 
ern edge of the maser distribution (Fig. H. The observed parame¬ 
ters of this continuum source imply a brightness temperature of 
~5400 K. For a linear maser geometry dWatson & Wvldl2003l) in 
unsaturated regime with R/T < 0.2, where R is the pump rate and 
F is the decay rate, a modulation of the background emission can 
produce the maser flux density variations (Sobolev & Parfenov, 
priv. comm.). The disappearance of cyclic changes and weak de¬ 
cay of the 97.10 km s“' fea ture after MJD 56208 resemble those 
reported for G188.95+0.89 dGoedhart et al.ll200^ which has been 
interpreted in terms of the recombination of the ionized gas against 
which the maser cloud is projected. 



Figure 7. Phase-folded normalized light curves of selected features of 
four targets; a) G75.76+0.34, b) G73.06+1.80, c) G45.473+0.134 and d) 
G25.411+0.105. Each panel is labelled by the velocity of feature. 


4.3 Size of maser region showing periodic variability 

G22.357+0.066, G25.411+0.105 and G45.473+0.134 have been 
observed at 6.7 GHz with milliarcsecond (mas) resolution 
dSartkiewicz et alj2009Ll20l^ . Their maser cloud maps are shown 
in Figurel^together with the position of the nearest infrared source 
for which the spectral density distribution is consistent with an 
HMYSO candidate. The groups of the maser clouds with clear and 
marginal periodic variability are marked. 

In G22.357+0.066, the clouds showing flares are located in an 
elongated area of size <156x29 mas (Fig. 0 that corresponds to a 
projected linear si ze of < 760 x 14 0 au for the near kinematic dis¬ 
tance of 4.86 kpc dReid et alj|2009h . The maser clouds without or 
with marginal variations of the intensity are in the northern edge of 
the distribution of size >1500 au. VLSI observations at two epochs 
spanning two years have shown that the spatial structure of th e 
source on a mas scale remains unchanged dSzvmczak et al.ll201 ll) . 
The maser spectra obtained with the VLBI and single dish at nearly 
the same epoch are very similar. This means that almost all maser 
clouds have compact cores and the source structure was full y recov¬ 
ered with a 6masxl2mas beam dSartkiewicz et alj|2009l) . There 
is a 5 d delay in the peak of the flares between the 80.09 and 
83.42 km s“* features (Table[T]l that can be explained by light travel 
time between the maser clouds of 80.23 and 83.57 km s“* (Fig. 0 
that yields a distance of 850 au between the regions. This agrees 
within the error range with the projected distance between these 
clouds of 700 au. Thus, the geometry of maser clouds and time de¬ 
lay measurements clearly indicate that an exciting source is located 
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Figure 8. Histogram of periods for periodic 6.7 GHz methanol maser 
sources. The objects with asymmetric flare profile and sinusoidal-like pro¬ 
file ai‘e shown in blue and red, respectively. The conti'ibution from the Torun 
32m telescope observations are marked by the red line. 
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Figure 9. Maps of the 6.7 GH z methanol maser in three periodic sources 
iBartkiewicz et alJl20Qgl[^l4) . The circle size is proportional to the loga¬ 
rithm of maser brightness of cloud and its colour corresponds to the veloc¬ 
ity. The individual maser clouds are also labelled by the velocity in km s“^. 
The groups of maser clouds showing clear and marginal periodic varia¬ 
tions are framed by the solid and dotted rectangles, respectively. The star 
symbol with error bars in each panel indicates the loca tion of the nearest 
infrai'ed source for G22.356+0.066 and G45.473+0.134 jLucas et al.l2Q08l) 
and G25.411+0.105 foe Buizer et al .120121) . The triang le symbol wife erro r 
bars shows the position of the radio continuum source iHofner et alj201 j]) . 


in the west side of the maser distribution, consistent with the lo¬ 
cation of the infrared HMSYO candidate marked in Figure|3 The 
brightness temperature (Tb) of the maser clouds in G22.357-l-0.066 
derived from Bartkiewicz et al.’s data ranges from 6 x 10^ to 
4 X 10* K. The mean values of Tb for the clouds with periodic and 
non-periodic emission are 1.7 ± 0.9 x 10* and 1.3 ± 0.3 x 10* K, re¬ 
spectively. A cluster of the brightest maser clouds showing periodic 
variations is located closer to the exciting source than the weaker 
clouds with or without periodic behaviour. 

The above summarised observational characteristics of 
G22.357-1-0.066 do not support the hypothesis that the maser flares 
are related to changes in the background free-free emission. The 
maser region with flaring behaviour has a size more than one order 
of magnitude smaller than the typical si ze (<10000 au) of hyper¬ 
compact Hit regions (iHoare et al.ll20()7h . No obvious mechanism 
exists to create such small scale fluctuations in the background ra¬ 
diation. Moreover, as the clouds showing periodic variations are 
brighter than those without or with marginal variations one can sup¬ 
pose that the ratio of RjT is higher in the first group of clouds than 
in the second group. Thus, the flux density of flaring components is 
expected to be less dependent on the variations of the background 
photon flux. 

The observations of G22.357-)-0.066 can support the hypothe¬ 
sis that the maser flares are related to changes in the maser optical 
depths and excitation temperature in some maser regions periodi¬ 
cally triggered by a pulse from the massive exciting star or binary 
system. Periodic pulses of UV radiati on are predicted in th e model 
of a pulsationally unstable protostar dlnavoshi et alJlTOlsh but it is 
difficult to explain why only some parts of the maser regions in 
G22.357-1-0.066 are excited. The same difficulty applies to the sce¬ 
nario where periodic variations of the dust temperatur e are caused 
by va riability in the accretion rate in a binary system dArava et al.l 
120101) . The illumination of the molecular disc by the bow shock 
hot and den se material of the spiral sho ck in the circumbinary disc 
central gap dParfenov & Soboi^l2014h could explain the periodic 
maser flare of some part of maser regions, in a particular condition 
when the maser emission arises in a disc seen edge-on. 

We conclude that any mechanism responsible for the periodic 
changes of the maser emission must be able to influence only part 
of the masing region if it is to be consistent with observations. 

VLBI observations of G45.473-)-0.134 in March 2010 have 
shown that the maser emission comes from three clusters 
dBartkiewicz et alj|2014l) . No significant changes in morphology 
were seen as compared w ith the maps obtained in 2007 with MER¬ 
LIN dPandian et alj20ll]) . Comparison of the spectra taken with the 
32 m telescope and VLBI indicates that less than 15 per cent of the 
flux d ensity was missed in the VLBI observation dBartkiewicz et al.l 
|2014 their figure Al), thus the source morphology is accurately 
recovered. The periodic flares originate from a single cloud at 
59.60 km s y‘ (Fig.l^of Tb = 7.5 x 10* K. In the map taken with 
MERLIN dPandianet al.ll201lh the maser distribution has a dou¬ 
ble structure of size less t han 70 mas that corr esponds to 480 au 
for the distance of 6.9 kpc dPandian et alj|20()^ . The cloud show¬ 
ing the periodic emission is separated by 413 mas (2850 au) from 
the cluster of four clouds at velocities of 64.29-66.30 km s“*. This 
cluster shows weak periodic emission only during some cycles well 
synchronised with the flares of 59.60km s“' emission. No periodic 
emission is seen from the 62.30 km s“* cloud that is located 1.03” 
(7100 au) south-west from the site of the periodic emission. We 
conclude that the mechanism causing periodic changes of the flux 
density is effective only in a relatively small (<480 au) region of 
the source. As the position uncertainty of the candidate of HMYSO 
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observed in the infrared (Fig. is comparable with the size of the 
overall maser distribution it is plausible that it may drive the maser 
periodicity. 

G25.411+0.105 shows a ring-like distribu tion of the maser 
cloud s (Fig. 1^ with a diameter of 190 mas llBartkiewicz et al.l 
llOOgl) that correspond s to 1000 au for the distance of 5.35 kpc 
jPe Buizer et al .l2012h . The maser emission does not show any ve¬ 
locity gradient; some blue- and red-shifted clouds are considerably 
blended spatially and in some spectral channels (e.g. 94.82 kms“') 
the emission appears at very different parts of the ring. In such a 
case the light curve of the blue-shifted emission can significantly 
differ from that of the red-shifted emission as we observed with 
the single dish. The range of brightness temperatures of the maser 
components is 0.3 - 12 x 10* K. The unresolved ( 0.8”x <0.6“) con¬ 
tinuum source at 25.5 GHz jHofner et alJbOlU) coincides within 
200 mas with the strongest maser component. This continuum 
emission is interpreted as an i onized jet at the ba se of a massive 
flow observed in the CO lines jHofner et alj|201 ih and it may in¬ 
duce the maser variations. If the source could be detected at cm 
wavelengths it would be possible to evaluate the effect of the back¬ 
ground emission on the maser flaring. 

We have inspected the first three objects listed in Table 
[J for continuum emission at 5 GHz using the Cornish survey 
i Hoare et alj[2012l; IPurcell et al.ll2013l) . No emission was seen at 
the 5cr level of 2mJy. We found no evidence for cm radio emis¬ 
sion in the remaining targets when searching the literature using 
the SIMBAD database. One can speculate that our periodic sources 
are associated with HMYSO in binary systems in very early phases 
with vigorous accretion processes. 


5 CONCLUSIONS 

We have found four new periodic methanol sources at 6.7 GHz 
and updated the light curve of the previously known source 
G22.357+0.066. This study increases the number of known peri¬ 
odic masers to 16 with a median period of 210 d. The period distri¬ 
bution hints to an association of the periodic masers with massive 
binary systems. The diversity of variability patterns cannot be eas¬ 
ily interpreted within the frame of a simple model of the colliding- 
wind binary system which provides periodic variations of the con¬ 
tinuum emission from the background Hii region. The study re¬ 
vealed that the periodic variations can be confined to a single fea¬ 
ture in the spectrum for some sources. In two sources, the maser 
flares originate in a region of linear size less than 480-760 au that 
contributes 20-50 per cent of the overall maser distribution. This 
implies that the mechanism generating periodic maser variations 
affects only some parts of the masing regions. 
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Figure Al. False-colour image of the 6.7 GHz methanol maser flux density Figure A2. Same as in Fig. lAll but for source G25.411-t-0.105. 

over velocity and time for G22.357+0.066. The velocity scale is relative to 
the local standard of rest. The horizontal bars in the left coordinate corre¬ 
spond to the dates of the observed spectra. 
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APPENDIX A: T HE DYNAMIC SPECTRA AND L-S 
PERIODOGRAMS 

For completeness, in this appendix, we present the dynamic spectra 
of the sources and the Lomb-Scargle (L-S) periodograms for the 
features for which significant frequencies were found. 
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Figure A6. The L-S periodograms for the features of the sources for significant frequencies, i.e. with a false alarm probability < 10“*^. The panels are labelled 
with the source name and the velocities of features considered are given. The dotted horizontal lines mark a false-alai'm probability of 0.01 per cent. Note 
harmonic series in the periodograms for the features of sources G22.357+0.066 and G45.473+0.134. These significant frequencies seem unlikely to be real 
detections. 
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Figure A3. Same as in Fig. lAll but for source G45.473+0.I34. 
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Figure A4. Same as in Fig. lAll but for source G73.06+L80. 
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Figure AS. Same as in Fig. lAll but for source G75.76+0.34. 




















